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Summary 


In this report, an analysis is presented of the acoustic field which 
is generated by a circular jet in a weak crossflow. The main portion of this 
report contains the analysis for the acoustic intensity which would be meas- 
ured from the jet flow in an ideal wind tunnel (i.e., a wind tunnel with no 
reflecting walls). For this situation, the movement of the acoustic medium 
results in the acoustic intensity being greater upstream of the jet flow than 
at ah equal distance downstream. The Appendix of this report contains a 
summary of the analysis of the acoustic intensity of a crossflow jet in 
motion such as would be generated by the lift jet of a V/STOL aircraft under 
take-off conditions. The result of this work is a comparison of the direc- 
tional distribution of acoustic intensity for an ideal wind tunnel test of 
a jet and a corresponding flight test. 
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SYMBOLS 


speed of sound in the acoustic medium 

r 2 L 2-,% 

Cx l + x 3 ] 

diameter of jet aperture 
function of separation distance 
acoustic intensity 

acoustic intensity per unit volume of turbulence 
characteristic length in the turbulence 

convection Mach number (a M is the corresponding velocity vector) 

o — 

V. 

— ^ — Jet Mach number 


— crossflow Mach number or aircraft Mach number (a Q N is 

o the corresponding velocity vector) 


P ijkl^’-’ x) stress correlation function 
P thermodynamic pressure 


P. . = 

ij 


3u, 3u . . 3u, 

PS. . + y C + 

xj 3x 3x 3 Sxj^ xj 


stress correlation function 


[(1 - N 2 )jx - yj 2 + N 2 (x- - y,) 2 ^ - N(x - y_) 


1 - N 


2 J _r 

.,2 


r a - N (x 2 - y 2 ) 
r g - N (x 2 - z 2 ) 

(y) Cy + h) second-order velocity correlation function 




iii 


ij 


' J 


x - z 


[(1 - N 2 )|x - z\ 2 + N 2 (x - y,) 2 ^ 


P ij + pu i u j " a o p<S ij 

time (in different coordinate systems - t, t) 

I* “ 2.1 

t 


o 

lx - z I 


t - 


U, 


u. 

X 


mean velocity component 


- velocity component in turbulent flow 




turbulent velocity component 


x 


U i X i 


- component of velocity in the direction of the observer 


V J 


u^ + ~ total velocity component 


z . 
x 


6 0 ^) 


jet exit velocity 
crossflow velocity 


spatial coordinate (in different coordinate systems - x , 


x ±J - x i ) 

position vector of the observer 

spatial coordinate in the turbulent flow 

position vector in the turbulent flow 

spatial coordinate in the turbulent flow 

position vector in the turbulent flow 
1 


d - N 2 ! 5 

Dirac Delta function (, 6 (x^) = o,x^o; j 6(x^)dx^ = 1, a<o<b) 

* a 


r b 


iv 



A 


e 

n 


0 

x 

_x 

y 

? 


a 

T 


T = 

a 

T e 

<t> 

ij 

“f 

Subscripts 


V 

\ 

(q - i £ j ~ t ^ ie Kronec ^ er Delta function 
z_ - 2 . - separation vector 

distance between elements of the quadrupole source 

position vector in the turbulent flow (in a moving coordinate 
system) 

polar angle between the observer and the x^-axis 

111 


A “ ~ separation vector 

dynamic viscosity 

position vector in the turbulent flow (in a moving coordinate 
system) 

total density 

equilibrium density 

parameter related to eddy size 

(t„ - t ) or (t„ - t ) - separation time 
6 a 6 a 



t - R /a - time of emission of signal 
a o 

t - R 0 /tf - time of emission of signal 
p o 

polar angle between the crossflow direction and the direction 
of the observer 


function related to p such that p 


3 2 <j>. . 
p = —il 

K o 3x.3x. 
i J 


characteristic frequency in the turbulence 


indices taking on values of 1,2,3 
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I. Introduction 


The acoustic performance of aircraft components has become an important 
consideration. Recently, the wind tunnel has been used as a test facility to 
evaluate the acoustic properties of aircraft components under simulated flight 
conditions (ref. 1). In particular, lift fans for use in V/STOL (vertical/ 
short take-off and landing) aircraft have been tested in the 12.3 meter by 
24.6 meter (40 feet by 80 feet) wind tunnel of the Ames Research Center of the 
National Aeronautics and Space Administration. The lift fan is placed so that 
it exhausts in a direction perpendicular to the direction of flow in the wind 
tunnel. This arrangement simulates the take-off or landing configuration of 
a V/STOL aircraft in flight. The acoustic measurements are made from micro- 
phones which are placed at various positions around the lift fan. 

While tests in a wind tunnel of aerodynamic properties can be used 
directly to infer the properties under flight conditions, the acoustic data 
gathered in wind tunnel tests do not completely simulate the acoustic con- 
ditions obtained in flight tests. One source of noise from the lift fan 
may be considered to be the turbulent exhaust flow. When the noise from this 
source is measured in the wind tunnel, both the source and the observer (i.e., 
the microphone) are at rest, but the acoustic medium (i.e., the flow in the 
wind tunnel) is in motion. Signals which are emitted by the source must 
propagate through a uniformly moving acoustic medium in order to reach the 
observer. On the other hand, the type of acoustic data which is normally 
wanted from flight tests is the noise level from an aircraft in flight which 
is measured by an observer at rest on the ground. In this case, signals 
which are emitted by the turbulent flow from the engines of a moving air- 
craft propagate to the observer through an acoustic medium at rest (assuming 
calm atmospheric conditions) . This is a different acoustic configuration 
than that of the wind tunnel test. 

In this report an analysis is presented of the sound field which is 
generated by the turbulent flow from a lift fan in a wind tunnel with no 
reflecting walls (i.e., an "ideal" wind tunnel). The flow from the lift fan 
is modeled by a cold jet flow from a circular aperture. The jet exit vel- 
ocity is subsonic, and the jet exhausts into a uniform crossflow. The 
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directional distribution (or directivity) of the acoustic intensity which is 
calculated is compared with the corresponding directivity for a flight test. 
The results provide a comparison between the acoustic data taken in a wind 
tunnel and the corresponding data taken in a flight test. 


Characteristics of a Circular Jet in a Crossflow 

The model for the actual flow from a V/STOL lift fan is the turbulent 
jet flow from a circular aperture. In order to simulate take-off conditions, 
the jet will be considered to exhaust into a uniformly moving crossflow with 
velocity in the positive x^- direction (fig. 1). The exit velocity from 

the jet, V_. , is subsonic, and the jet flow is considered to have the same 
density as the crossflow. The resulting flow differs in several respects from 
that of the circular jet with no crossflow. There are three differences which 
influence the acoustic analysis . 


The most obvious feature of the crossflow jet is the deflection of the 
turbulent flow in the direction of the crossflow. The severity of the deflec- 
tion is related to the ratio of the jet velocity to the crossflow velocity. 

The position of the centerline of the jet, as determined by measurements of 
the maximum velocity in the jet, is indicated on fig. 2 for several values 
of ( re f* 2). It should be noticed that for V j /^2 as low as 8, the 

centerline of the jet is only deflected two jet diameters in the first eight 
jet diameters downstream from the aperture. In the following, a crossflow 
for which V j /^2 ^8 is referred to as a weak crossflow. 


Another prominent feature of the crossflow jet is the distortion of 
the initial circular cross-section of the jet as the flow proceeds downstream. 
The shear produced by the crossflow initially distorts the circular flow into 
one with a cross-section of "kidney" shape (see fig. 1). As the flow proceeds 
further downstream, a typical cross-section contains two counter-rotating 
vortex-like structures (see ref. 3). 


Finally, in the immediate vicinity of the jet flow, the crossflow is 
deflected by the jet. In this regard, the jet flow has been likened to a 
porous cylinder which acts as an obstacle to the crossflow (see ref. 3). 


2 


Acoustic Analysis 

The conventional use of a wind tunnel is to conduct tests which simu- 
late flight conditions. In order to do this, the air in the wind tunnel is 
moved uniformly past a test model. When acoustic measurements are made in a 
wind tunnel, signals which are emitted from a source must propagate through 
the moving medium to reach an observer (e.g., a microphone). Acoustic tests 
performed outdoors, however, generally involve measurements made on the 
ground of the noise generated by the engines of an aircraft either in flight 
or on the ground. A fundamental difference between acoustic measurements 
made in a wind tunnel and outdoors is that in the wind tunnel the acoustic 
medium is in motion relative to the observer while outdoors the medium can 
normally be considered to be at rest. 


The analysis of acoustic propagation in a moving medium begins from a 

consideration of the continuity and momentum equations. The acoustic medium 

is assumed to be uniformly moving with a velocity of in t ^ ie positive 

x^- direction (see fig. 1). The total velocity at any point in the 

medium is the resultant of the convection velocity and other velocities 

u^ which might result from an acoustic disturbance, turbulence, etc.; thus, 

V. = u. + V„ 5.„ where is the Kronecker Delta (6 = { S* ^ ). 

1 1 2 12 x2 x2 1 x = 2 

The continuity equation is written 


an 

le. + = o 

at 3x. 


or. 


3p 3pu i. 3p 

-Til 4. — + V -rr_ = n 

at 3x. 2 3x_ 

x 2 


( 1 ) 


since is uniform. 


The momentum equation is written 


*PV ± 

at 


3p7.7. 3P . . 

z_j_ hi 

3x . 3x . 

3 3 


where, for a Newtonian fluid. 


P. . 
13 


. 37. 

= P6. . + p C- - 3 — 

xi 3x. 

J 3 


^j2 

3x^ 3 3x^ 


6 ..] 

13 
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Upon substitution for V and use of the continuity equation, the momentum 
equation is obtained in the following form: 


3pu ± 

3t 


3pu.u. 

U 

9x . 

3 


Spu . 

+ V — — = 
2 3x. 


3P . . 
.. U 
3x . 

3 


( 2 ) 


The continuity and momentum equations are combined by taking the partial 
derivative with respect to time of eq. 1 and subtracting from it the deriv- 
ative with respect of x. of eq. 2. Use of the continuity equation and 

1 2 
2 3d 

subtraction of the term a r from both sides of the resulting equation 

0 3x. 

l 

gives the following governing equation: 


2 2 

UL+ 2V O— 

at 2 2 3t3x 2 


+ v, 


2 9 2 P 

1 3 2 

X 2 


2 3 P _ 


3x. 

l 


3 2 T. . 

il 

3x. 3x. 
i 3 


(3) 


where T,, = P . . + pu.u. - a p6,. . The left-hand side of eq . 3 describes 
ij ij i 3 o p ij 

the propagation of acoustic waves in the uniformly moving medium. The right- 
hand side is a quadrupole source involving the stress T _ which results 
from pressure, viscous stresses, and momentum flux. In classical acoustics 
this stress is negligible. In a region of turbulence, however, T is non- 
negligible primarily because of the momentum flux contribution. 


The solution of eq. 3 is obtained by using a method of coordinate 
transformations. This technique consists of finding a system of coordinates 
in which both the medium and the source are at rest with respect to the 


observer. In order to begin, a new function 


to the density as follows : 


ij 


" P o = 3^ 


is defined which is related 


(4) 


where 


is the undisturbed density of the acoustic medium. When eq. 4 


is substituted into eq. 3, the following equation for <j> . . is obtained: 


3 2 <f>. . 

^ + 2V 


3t 


3 2 4>. . 

2 3t3x^ + V 2 


ij 

2 2 

2 „ 2 x 3 = 

2 a o n 2 


3x r 


3x . 

l 


T. . 
ij 


(5) 
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In general, the stress T is a function of position in a region of 
turbulence, and of time, t. The stress can be expressed using the 

Dirac Delta function as T„(t)6(x^ - y^)6(x2 - ~ where 

6(x^ - y^)6(x 2 - y2)<$(xg - y^) is equal to zero at all values of position 
(Xf, X2, x 3 ) except at the location of the stress (y^, y2, y^) . A solution 
to eq. 5 is first obtained for a stress located at the position (y^, y^, y^) . 


The equation to be solved is the following: 


3 2 4> . . 

+ 2V 

at 


s2 hj 2 
— ±J- + v. 


2 at ax, 


a 2 <j) . . 

n _ 


2 a 2 
8X2 


a 


ifhi . 

o a 2 
3x. 

1 


T_(t) 6(x 1 - y 1 )6(x 2 - y 2 )<S(x 3 - y 3 ) 


( 6 ) 


The left-hand side of this equation can be transformed into the normal wave 
operator on <(>„ by the following Galilean coordinate transformation: 


X 1 = X 1 ~ y l ’ x 2 = x 2 “ y 2 _ V 2 C 


x 3 = x 3 - y 3 ; t 


= t. 


The relationship between derivatives is the following: 


_a 

3x 


r^’ 


a 

3x„ 


a . _a_ 

ax 2 ; ax 3 


a a_ _ a a 

3 x 3 ; at at “ 2 a£^ 


The above relationship between coordinates is used to transform eq. 6 into 
the following equation: 

a 2 <j) a 2 <j> 

P- - a 2 ^4- = T (t)6(x 1 )6(x 7 +V„t)6(x_) . (7) 

3t 2 0 3x 2 12 2 3 

1 

This transformation has a simple interpretation. In the new coordinate 
system moving with velocity in the positive X2~direction, an observer 

finds the acoustic medium to be stationary. Relative to the observer, 
however, the source moves with a velocity in the negative X2~direction. 

The series of coordinate transformations by which eq. 7 may be solved 
is outlined in Morse and Ingard (ref. 4). The first of these, a Lorentz- 
type of coordinate transformation, has the property of leaving the left-hand 
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side of eq. 7 in the same form; however, the source term on the right-hand 
side becomes stationary. This coordinate transformation is given as follows 


x, = x 


L ; ^2 = y ^2 + V 2^ ’ ^3 = *3 > ^ = Y(t + ^- x 2 ) 

o 


where 


Thus, 


3^ 39* 


y = (1-N 2 ) ^ and N = V 0 /a 

Z O 


3 , 3 , N 3 . 3 3 3 ,3 , „ 3 N 

3C = y( 357 + 7" 7? ); 357 = 357 ; lf = Y % + V 2 W/ 


' 2 2 ~o *'•' “~3 ““3 

In terms of the new coordinates, eq. 7 becomes 

2 2 

3 <t> . . , 3 

1 1 - a r 1 = yT. . (y^)S(x )S(x„)6(x ) 

Z O Z XJ i x J 


3 r 


( 8 ) 


3x. 

x 


The above form of the equation results from use of the property of the Dirac 
Delta function that d(x 2 /y) = y<$(x 2 ) • 

One final coordinate transformation is needed to put the governing 
equation into a form for which the solution is immediate. This transforma- 
tion is the following: 

'X 


f\j 'Xj ^ f\i 

1 = yx^ ; x 2 = y x 2 ; x 3 = y x 3 ; t = yt . 


In terms of these coordinates, eq. 8 becomes 

,2 

2 ‘ 
a 


3 2 * 3“<(, 

' a o = y2 T (t)6(x 1 )5($ 2 )6(x 3 ) 


3t 


3x . 

l 


(9) 


This equation represents wave propagation in three dimensions through a 
stationary uniform acoustic medium resulting from a stationary source of 
strength ’ 
following: 


2 -*■ 

strength y T„ (t) located at the origin. The solution to eq. 9 is the 






( 10 ) 


4 m r 
o 


6 


where 


r = 


& 


+ X„ 


-*■ 2-, 

x 3 J 


It should be noted that in order to determine the function 


position r and time 


t, 


3-3 


at any 


the acoustic signal, t - — 


13 


must be evaluated at the time of emission of 


a 


The independent variables (x^, t) in the final system of coordinates 
are related through the coordinate transformations to the independent vari- 
ables (x^, y^, t) in the original coordinate system; thus 


r 


[(1-N 2 ) |x-yj 2 + N 2 (x 2 -y 2 ) 2 lT 2 


1-N 


and 


where 



R = 


C(l-N 2 ) |x-yj 2 + N 2 (x 2 -y 2 ) 2 J^ - N(x 2 ~y 2 ) 
1-N 2 


(ID 


In terms of the original coordinates the solution for 



is the following: 


3-3 


(x.t) 


T.,(t - ?-) 


4tt a Q C(l-N ) |x-yj + 



( 12 ) 


This equation determines the value of <)>.. which results from a single 
source T„ located at the position (y^, Y 2 > Yg) • When a region of turbu- 
lence is present, the solution for <j>„ is generalized by integrating the 
result from each source in the region. Thus, for a distributed source. 


cj>_ (x,t) = 




4ir a Q 2 C(l-N 2 ) |x~x| ^ + N 2 (x 2 -y 2 ) 2 H 2 


d 3 x 


(13) 


where is the position in the region of turbulent flow, and the integral 
extends over the entire region of turbulence. 
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The acoustic density disturbance at a position x and time t is 
related to 4>_ by eq. 4 and is therefore determined by the following 


, .x 13 

P(-’t) - P 0 = : — 2 3^177 

4tt a x j 

o 


T ±j <z>t - f- ) 

j 1 % .3 

o o o o U ® Z." 

C(l-N ) |x-2_| W(x -y rj* 


In this equation the retarded time (i.e. the time at which the signal is 
emitted) is t - — where R is defined in eq. 11. The quantity R 
represents the distance that the signal propagates with speed a through 
the moving medium from the source to the observer. This interpretation can 
be demonstrated with the aid of fig. 3. The plane indicated in fig. 3 
contains both the observer located at position (x^, x^) and the source 

located at the origin. The direction of the observer from the direction of 
the crossflow is given by angle <j> . In order for a signal to be received 

by an observer, it must propagate through distance R = At where At is 
the time interval between emission and reception of the signal. During the 
time interval At , the signal is convected a distance = NR by the 

acoustic medium. These relationships allow R to be evaluated in terms 
of r = |x] and <|> as follows: 


R^ = (rsin<j>)^ + (rcos<j> - NR)^ , 

and therefore, 

R = 

When N is less than unity, the positive sign is taken, and the relationship 
rcos(|> = x^ is used to obtain the value of R in eq. 11 for a source located 
at the origin. 



Equation 14 provides the formal solution for the acoustic density dis- 
turbance when the acoustic medium is moving uniformly with velocity in 

the positive x„-direction. When the Mach number of the medium is such that 
2 1 

N «1, the denominator of eq. 14 is approximately equal to the distance 


between the source and the observer. In this limit, however, the retarded 

I -~¥ I 

txme is not given by t - J *- as is the case for a stationary acoustic 
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medium. At small Mach numbers, the primary effect of the movement of the 
acoustic medium is therefore the alteration of the time between emission and 
reception of acoustic signals. The net acoustic signal which is received at 
a point from sources located near to each other is strongly dependent upon 
the time interval between emission and reception of the signal from each 
source. The alteration in the retarded time by the movement of the acoustic 
medium therefore influences the acoustic field generated by nearby sources 
(e.g., the acoustic field which is generated by a region of turbulent flow). 
This effect can be put in a quantitative form by carrying out the differ- 
entiation of the integrand of eq. 14 as follows: 

D , 

‘ . 0. . 1.. J. j. (JA, U, Vfc - . 

J * * * J O J 


V _ 


■) = 




3r. n 3T. . 


3R 


where 


r A = C(1 - N 2 )|x-yJ 2 + N 2 (x 2 -y 2 ) 2 H' s 


and 



Therefore, 

a 2 

3x.3x. r. 

i 3 * 


1 3R 3R 

2 3x. 3x. 
a r. x i 

o * J 


3 2 R 


3T 


3x, 3x . ' 
i 3 


il 


3x ’ 


3 T ij _1 , 1 _ 3R_ 1_ _3R_ 

a T ’ 2 V* r * 8x i 8x j r * 3x ± 3x j 

1_ 2 _^*_^±_ 92r * 

+ r * ^ r * 3 x ■ 3x - 3x. 3x . ^ T ij 

**13 i j J 


When the solution is restricted to the acoustic far field (i.e., to distances 

from the source which are large compared to (2 tt) ^ times a typical acoustic 

wavelength) , the dominant term 0 of the integrand is the one that decays most 
^ * * —1 
- 1 — 1 .. tl. — -r 11 — j — - while the coefficients 

-2 , _ -3 


.,2 


is of order r. 


slowly. The coefficient of 

3T 1 . 3t 

of r i " and T . . can be shown to be of order r 
3t 13 


and r . 


respectively; 


therefore. 


3x. 3x. 
1 3 


T. . (£ t’) , 

, ij ^ . _ 1 3R 

v _ ) o a-u- 


2 

a r A 
o * 


3R • 
3x . 3x . 


3 2 T. 


3t 


11 

,2 ' 
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The derivatives of R are easily calculated with the result that 


3R 

3x 


1 (1-N j 


i (x-y) (1-N 2 ) + (x -y )N 2 6 

1 [ — i — i ±_± — _ N6 ] 

2, u r x i2 J 


The density disturbance in the acoustic far field is therefore given as follows: 

2 2 


p(r *> t >- p o ‘ 4~ X 2 

4ira (1-N ) 
o 


r 2 ^x.-y.Ml-N ) + (x i -y i )N 6 ±2 


- N6 i 2 3 


(15) 


(x-y )( 1-N ) + (x -y )6 . 3 T. (y_» T ') 3 

[— ■— 1 . ■ J J j 2 - *«.,] — ^4 


J2- 


3t 


,2 


A somewhat unintuitive result is contained in eq. 15. When the source 
is a quadrupole with either axis in the direction of the crossflow (i.e., 
with either i or j equal to 2), a greater density disturbance is calculated 
upstream of the source than at an equal distance downstream for all observer 
positions except along the x 2 -axis. This result can be demonstrated by con- 
sidering the source to be a single lateral quadrupole T^ 2 located at the 
origin. This type of source represents a shear-stress on the x^-plane in 

the x„-direction or vice versa. For this example the square of the crossflow 

^ 2 
Mach number is considered to be small compared to unity (i.e., N «1) . The 

density disturbance is given by eq. 15 to be the following: 


P(r,t) - P Q = 


4it a r 
o 


X X 

4~ <!T )( r " N) 


2 

3 T 


3t 


12 _ 

,2 


, where r = |x| 


(16) 


A stress is thought of as resulting from two equal and opposite forces acting 
on a fluid element. Each force can be represented by a dipole (i.e., an 
adjacent source and sink of equal strengths which emit simultaneously) . There- 
fore, the quadrupole source T^ 2 is represented as shown on fig. 4 where the 

+ indicates a source and the - indicates a sink of equal strength. In the 
absence of a crossflow (i.e., with N = 0), the directional distribution of 
the density disturbance at a fixed distance r which is generated by the 
X 1 X 2 

source is (^— ) (— ) . This pattern is the symmetric "rose" shown on fig. 4a. 
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The reason that an observer along the x^- and axes finds zero density 
disturbance is the interference of the signals which are simultaneously 
emitted from the sources and sinks comprising the quadrupole. The observer 
on either axis is equidistant from source-sink pairs of the same strength; 
thus, the signals which are simultaneously emitted from each source-sink 
pair of the quadrupole propagate through the medium with velocity a and 
arrive at the observer position simultaneously. The result of a signal from 
a source arriving in phase with a signal from a corresponding sink is complete 
destructive interference and zero net density disturbance. For the observer 
located off-axis, however, the corresponding source-sink elements of the 
quadrupole are not equidistant from the observer and the interference is not 
complete. When the acoustic medium is moving with a velocity < 2 q N, signals 
which propagate in any direction are simultaneously convected with the medium 
as they propagate. Convection by the acoustic medium has the effect of alter- 
ing the time required for the signals from each element of the quadrupole to 
reach the observer. This gives rise to a distorted distribution of the acoustic 
density disturbance. This directional dependence is given by eq. 16 to be 


x x 

(— ) (— - N) . 

r r 

Figure 4b gives the directional dependence of the density disturbance for the 
lateral quadrupole T ^ which results when the medium is moving in the posi- 
tive X 2 ~direction with a Mach number of N = 0.25. Two features distinquish 
this pattern from the symmetric pattern of fig. 4a. First, zero net signal 
is found not along the x^-axis but rather along a ray downstream of the 
x^-axis. This is explained directly by considering the effect of convection 
of the medium on the propagating signals . The signals which would have 
destructively interfered with each other along the x^-axis when there was no 
movement of the medium are now convected downstream. The second feature of 
the distorted "rose" pattern is the enhanced signal received upstream of the 
source compared to that received at corresponding positions downstream. This 
feature can be obtained by considering the propagation of signals from each 
element of the quadrupole to an observer located at an angle 6 from the 


x-^-axis (measured counter-clockwise) and distance r from the origin. The 
separation distance between adjacent sources and sinks in the quadrupole is 
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e is small such that (—)<<!. Each element of the quadrupole (i.e. 


from the observer, where 


e , where 

each source or sink) is then located a distance 

the index i refers to the quadrant on which the element is located (e.g., the 

source on the first quadrant is located a distance SL^ from the observer at 

position (r,0). The observer distance from the source elements of the quad- 

- e 


^ 3 = r [l + — (sin6 + cos0) ] 


rupole can be calculated to be the following: 

Similarly, the observer distances from the sink elements of the quadrupole are 
^ 4 = r ^ ^ 2r (sin0 - cos0)]. Since the velocity of the medium is a Q N on 
the x 2 ~direction, the resultant propagation velocity in any direction 0 is 
a (1 + N sin0) . The time required for a signal to traverse a distance H. 

° &i 1 

is therefore given approximately by t^ = — r^y • (Note, the approx- 
imation is that the signal from each element propagates on the direction 0.) 
The source elements of the quadrupole may be considered to fluctuate with a 
frequency o> as sin rnt while the sink elements fluctuate as -sin wt. A 
signal which is emitted by any element of the quadrupole at time t arrives 
at the observer at a later time t + t^; for instance, the signal which is 
received by the observer from the sink located on the second quadrant is 
- sin [u)(t + t 2 )]. The net signal received (NSR) by the observer is then 
calculated by summing the contributions from each element of the quadrupole 
as follows : 


NSR = sin [w(t + t^)] ~ s ^" n ^(t + t 2 ) ] 


This calculation for 
the result that 


+ sin [w(t + t^)] - sin [w(t + t^)] • 

NSR is simplified by using the smallness of 


NSR = 


sin9 cos0 


,u)e.2 . r 

(—) sm Lw(t + 


< 7 > 


with 


a 


The amplitude of NSR 


(1 + N sin0 ) ‘ 
is given by the factor 


a (1 + N sin0) 
o 


)] • 


sin9 cos0 ( 0 )E) 2 
(1 + N sin0) 2 a o 


which 


provides the approximate directional properties of the quadrupole source. The 

symmetric "rose" pattern (i.e., sin0 cos0) is modified by the convection 

_2 

factor (1 + N sin0) to provide an enhanced upstream NSR compared with 
that at corresponding downstream positions. This modification occurs because 
the destructive interference of the signals from each element of the quadrupole 
source is less complete for observer positions upstream of the source owing to 
the reduced propagation velocity in that direction. As a result, the upstream 
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"petals" of the distorted "rose" pattern are larger than the downstream "petals". 
It should be noted that the greater effectiveness of quadrupole sources upstream 
than downstream also applies to simple sources (monopoles) and to dipoles when 
the axis of the dipole is aligned with the -direction. 

Human response to a sound field is determined more by the acoustic 
intensity than by the acoustic density disturbance. In the acoustic far field, 
the intensity is defined to be the average rate of flow of acoustic energy 
through a unit area normal to the direction of wave propagation. When the 
intensity is integrated over a large surface which encloses an acoustic source, 
the total power output of the source is obtained. The relationship between the 
acoustic density distruban ce and the acoustic intensity is given as follows: 

3 

1 = (P(x»t)-P D ) 2 

o 


where the overbar indicates the time-average of the quantity underneath. (This 

is an approximate relationship for the acoustic intensity for small crossflow 

Mach numbers . It is a reasonable form for the example calculations which are 

included in this report. The exact expression for the intensity is given in 

ref. 14.) The acoustic intensity in the far field is determined by substituting 

eq. 15 into the above relationship. In order to simplify the results that 

follow, it will be assumed that the square of the crossflow Mach number is small 

2 

compared to unity (i.e., N <<1) . With this assumption, the acoustic intensity 
is given as follows: 


I = 


2 5 

16ir a. p 

o o 


1 r V y i 


r a r 3 


r - N6 i2 ] ' N6 j2 ] 

a a 


x, -z, 

i Js zr L " N<5 k2 3 [ 


e 


x„ -z 


1 L - N6 _] 




p 2 J 2 

9t 

a 


9x 


k£ — 6 ,3 ,3 
d £d z_ 


3 


(17) 
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where y_ and ?l are the positions in the turbulent flow, r = |x ~ zl » 

R R B a ■ 

r B = I* - *|. T a = t - -f , Tg = t - — , R a = r a - N(x 2 -y 2 ), and R g = r g 

- N(x 2 -z 2 ). ° ° 

The integrand of eq. 17 contains the product of direction factors and a 
time-averaged quantity which depends on the properties of the turbulence at 
different points in the flow field. It is useful to recast the integrand into 
a form which recognizes the eddy structure of the turbulence. This procedure 
is developed by Ffowcs Williams (ref. 6). First, a separation time t is 
defined such that t = x - x . The turbulence is specified to be a station- 

p ct 

ary function of time which allows the following relationship to be proven: 


, 2 4 R <18) 

3t i 3t B 3t 

* - . 

where the correlation function R - T ij ( ^ T „ )T k e <2.- T „ +T >- 

The above correlation function is determined by taking the time-average of the 

stress in the turbulence at two positions £ and and at corresponding times 

t and t + t . It is expected that the correlation function will have a 
a a 

maximum for positions ^ and jz which are close together and zero for positions 
which are widely separated. The region of high correlation corresponds to a 
turbulent eddy. 


The independent variable z^ is eliminated in favor of a separation 

3 3 

vector A^ which is related to £ and z_ as follows: A^ = z_-y_; d z_ = d i. 

k 

A new correlation function is defined by expressing R . in terms 

1J KJ6 

k 

of the variables A_, and t ; thus, R ijk£^- , “ ,T ^ = R ijk£^*“ ,T ‘ The 

separation time x is dependent on the observer position x and the positions 

R a -R B 

in the turbulence. This relationship is x = x„-x = . The distance 

B o a 0 

R a is a function of x and £ while R^ is a function of x and z_. The 

definition of the separation vector A allows R 0 to be expressed as a 

~ p 

function of x, y_, and A . Since the correlation function R — ^CZ’A* 1 ) 
is expected to be non-negligible only for positions in the turbulence which 
are relatively close to each other, the expression for x need only be 
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determined for values of A which are small compared to ( x - y ) . Thus , for 
an observer far from the turbulent field, there is only a small difference in 
(x - y) and (x - z) over the dimension of a turbulent eddy, and these dif- 
ferences can be neglected in the calculation of the acoustic intensity. The 

A.. (x-y)_NA 2 

expression for t is then t = ■ = — i — i . This expression gives the time 

a o ' — a o 

difference between the emission of signals from points located a small distance 
A apart. 


In a turbulent flow there is, in general, a velocity associated with the 
convection of a turbulent eddy in the flow. For instance, in the mixing region 
of a cold subsonic jet, the convection velocity of a turbulent eddy is approx- 
imately 0.65 times the jet exit velocity (see ref. 7, pp. 355-357). In order 
to estimate accurately the acoustic intensity, it is important to minimize the 
changes in R. .. 0 (£>A> T ) with respect to r which result merely from eddy 
convection. This is possible by defining a new separation vector X_ in a 
coordinate system which moves with the convection velocity a^M. Th e relation- 
ship between A and X_ is A = A + <3 0 H T • The separation time x is then 
determined as a function of X_ . In order to facilitate this transformation 
and the following results, it will be assumed that the convection velocity is 
of magnitude a Q M^ in the positive x^-direction only. This corresponds to 
the assumption of a weak crossflow. The resulting expression for x is then 

A- (x-z)- ' "2.1 


x = 


a o t|i“Zl " M l (x l“ y l )} 


A new correlation function P (£>A» T ) I s defined by expressing A 
terms of X and x ; therefore, T ) = R ijk£^-’A ,T ^ • Finally, the 

relationships between X_, A> and x are used to calculate the changes with 
respect to x of R^^^C^A.x) and P^j] c ^(2.»A* T ) • As shown by Ffowcs Williams 
(ref. 6) these two derivatives are not simply related because of the dependence 
of x on A and A • Th e relationship can be shown to be given by the 
following: 


It • { t 



|2-Y 3 

•Ux-z 

-M 1 (x 1 -y 1 )D 8x 


- Vl 
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When integrated over A_, the second term on the right-hand side can be shown 

to give a negligible contribution, and it is therefore neglected. The differ- 
3 3 

ential volumes d A_ and d ^ are related by the Jacobian of the coordinate 
transformation as follows : 

3 -j (x, -y -, ) 

d 3 A = d\ Cl - - 1 , 1 I 1 ] 

— x-y_ 


-1 


The acoustic intensity (eq. 17) can now be expressed in terms of the 
correlation function P^ ^ (z.»A» T ) as follows: 


I = 


. 2 5 

16ir a. p 
o o 


r 2 [1 - 
a 


i v. -y . x . -y . 

itt —c^-^ - N6.,: c-* 1 — 1 - N6.„: 

r a 12 r a J 2 


IJll 


-] 


(19) 


c-V- - '"W 7T mi 


3t 


This integral is simplified by restricting the observer position x to be 

far from the turbulent field compared to the dimensions of the flow field. 

Under this limitation, |x - zl = l^.l = r and x^-y^sx^, and direction factors 

become independent of the integration. Equation 19 is made consistent with 

2 

the assumption that N << 1 by expanding the four direction factors appear- 
ing in the equation to give the following terms : 


£ - Cr - BS k2 ] ^ - 

X.X.X^X^ 

- JL ^4 C6 i2 x j *k*£ + X i 5 j2 x k X JJ + x i x j 6 k2 X * + X i X j X k 6 ^2 :i 


( 20 ) 


The final expression for the acoustic intensity is in the following very 
useful form: 


I(r) = 


25 2 3 

16 it a p r [1 - ~±~\ 
o o r 


( X. X . X. x_ 


r i\ X SL N |-f . r 

f 4 3 [i i 2 s jV< “iVVl 


+ x l x j i k 2 x » + ViVu 3 ’ ^4 P ijkl < *"i' t) 


( 21 ) 
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Mx -5 

The effect of eddy convection is given by the Doppler factor Cl H 

while the influence of the medium velocity a Q N appears in the numerator 
of eq. 21. 


Even with the assumptions made to date to simplify the calculation, 
the solution depends upon the knowledge of the correlation function 
^ijki^*— ,T ^ at eac ^ P oint tll e turbulence. This implies knowing the 
stress at one spatial point and one time in the turbulence and at all other 
points at the proper separation time. In order to proceed beyond this point, 
it is necessary to use very simplified models of the jet turbulence. 


MODELING THE TURBULENCE 

In principle, the acoustic intensity can be calculated from eq. 21 if 
the correlation function (y^ A_, t) is known in the turbulent jet flow. 

At best, this function is very difficult to determine either experimentally 
or analytically for any turbulent flow. As a result, a very simplified model 
for the correlation function is used to calculate the qualitative aspects of 
the acoustic intensity which is generated by the crossflow jet. The method 
which is followed is given by Ribner (ref. 8) who calculates the acoustic 
intensity which results from a turbulent jet with no crossflow. 

For turbulent flows of low Mach number, the stress T„ is approx- 
imated as follows: T. .ap u.u. where p is the density of the surround- 

ij o x J, o 

ing acoustic medium. This approximation is reasonable for those flows in 
which heating or cooling is caused only by compression or expansion of the 
flow (see ref. 5). In this approximation, viscous effects have been 
neglected; the only contribution of viscosity would be a slight damping of 
the acoustic field outside the turbulent flow and a small contribution to 
the stress inside the turbulent flow. The correlation function in eq. 21 
is then given by the following: 

P ijk4 (Z.A.T) = p o 2 u. ( Z , T a ) u . ( X , T a ) u k (2+ A, T a +T ) u £ (y+A , T a +T ) . (22) 

This correlation function is substituted into eq. 21 to give the following 
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equation for the acoustic intensity: 


I(r) = 


2 52 M 1 X 1. 5 3t 4 

16^ a r (1 ) 

o r 


[u u u 'u ' - 2N {u^u u 'u * + 
xxxx 2xxx 


(23) 


where u = - 


u.x. 
l l 


u u u 'u ' }]d\ d^y , 
x x 2 x — 

is the component of fluid velocity in the direction of the 


observer, u^ is the velocity in the positive x^-direction and the prime 
indicates the velocity evaluated at + A_, + x) while the unprimed 


velocities indicate evaluation at (y,x ) 

a 


The velocity u represents the 

total velocity in the turbulent flow and can therefore be expressed as the 

sum of the mean velocity and the fluctuating component as follows: 

where U x is the mean velocity and u^ is the fluctuating 

This and similar expressions 
x 

it i .... i . . . i 

u 0 


u = U + u 

XXX 


component which is defined such that $ = 0. 

x 


for u' 


.>2 and u' 2 


are substituted into the terms such as 


u u u 'u ' 

xxx x 


which appear in eq. 23, and the indicated multiplication is performed. 
Several terms in the resulting expression can be neglected. Since the cor- 
relation function is to be operated on by the derivative with respect to x, 
those terms which contain no dependence on x are disregarded. Also, terms 


involving triple velocity correlations (e.g., u u u ' ) are neglected 

X X x 

because they are either zero or small depending on how the turbulence is 
modeled (see ref. 8). The result of these simplifications is given as 
follows : 


uuu , u'=4UU , uu' 
xxxx xx xx 


xxx 


When the other terms in eq. 23 are similarly considered, the integrand 
becomes the following: 


uuu'u'+4UU'uu'- 2N!I2(U„U ' + U U ') u u ' 
xxxx xx xx 2x x 2 x x 


+ 2U U ' (u„u ' + u u 0 ') + u„u u 'u ' + u u u 'tL' H 
xx 2x x2 2xxx xxx2 


(24) 
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At this point, the assumption is introduced that the turbulence in the 
jet is isotropic. For this assumption to be strictly valid, the intensity of 

the turbulence (i.e., Hu D ) at every point in the flow must be the same in 
all directions. This state of affairs is not found in the circular jet with 
no crossflow (see ref. 9), and there is no reason to believe that the assump- 
tion will be strictly valid for the circular jet in the presence of a cross- 
flow. The assumption however, does allow the calculation of acoustic inten- 
sity to be greatly facilitated, and it is expected that the result will pro- 
vide a reasonable first approximation. 

In isotropic turbulence, it can be shown (see ref. 10) that the second- 

order two-point correlations of the velocity in the same direction 

(e.g., u^u^' ) are non-zero while the correlations of a velocity at one 

point with a velocity directed perpendicularly at another point (e.g., u^u^') 

are equal to zero. Therefore, when the direction of the observer is the 

x„-direction, the terms u„u ' and u u„ ' are equal to u u When the 

observer is in the x^- or x^-directions , however, the terms u 2^ x ' ancl 

u u ’ are identically equal to zero. Thus, for any observer position, 

X z 

~ A | A A f A A f 

U 2 U x " U x U 2 = U x U x C0S< f> 

where <j> is the polar angle between the crossflow direction (x2~axis) and 
the observer direction (see fig. 5). Substitution of these relationships 
into eq. 24 gives the following expression: 


a i u 'u ' + 41 u 1 i i ' - 2N[2(U u ' + U U ' + 2U U ' cos*) u u ' 

XXXX X Y X X / V X 7 XX ” XX 


XX XX 


2 X 


x 2 


X X 


X X 


(25) 


■ A A I /V 1 , A A A lAl 

+ u„u uu+uuu„u J. 
2xxx xx2 x 


The expression given by eq. 25 contains terms which are second-order 
velocity correlations and fourth-order velocity correlations. While much 
is known about the behavior of the second-order velocity correlation in 
isotropic turbulence, little information is avilable concerning the fourth- 
order correlation. As a result, the assumption is made that the joint- 
probability distribution of u and u ' is normal. This distribution gives 

X X 
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the probability of finding a particular velocity at one point along with 

a particular velocity u ’ at another point. While this assumption is known 

not to be strictly accurate (see ref. 11) it does allow a relationship to be 

made between the fourth-order velocity correlation u u u 'u ' and the 

x x x x 

second-order correlation u x u x ' • This relationship is the following (see 
ref. 11): 


A y\ a * a a 

VjV* = 


^ "" A -1 I A I A A f . A A 

Vj ■ Vl + Vk • U jV + V; 


,<* t 

X 


U A 


The fourth-order velocity correlations in eq. 25 are written, 
relationship , as follows : 


u u u u 
xxx x 


’ = 2(G u ’) 

XX 


using this 


and 


A A A . A . A A A | 

u„u u u = u u u« 
2 x x x x x 2 


u 


= 2(u u ') cos<j> 
xx r 


(26) 


The terms which have no dependence on t have been neglected in the 

above relationships. The second-order velocity correlation u u ' is a 

x x 

function not only of the separation distance A_ but also of the separation 
time t. It is reasonable to expect that the dependence on x is exponen- 
tial so that the correlation can be written as u x u x ' (A_,x) = R^OOexp (-u f t) (27) 
where R (A) = u (y) u (y+A) is the simultaneous second-order velocity 
correlation, and is a characteristic frequency of the turbulence. Then 

eqs. 23, 26 and 27 combine to give the following equation for acoustic inten- 
sity: , 


I(r) = 


P o m f 


2 5 2. 

if a Q r (1 — ) 


{2R xx (1 - 4Ncos<j>) + 


(28) 


f*xx [U x U x’ (1 ” 2Nc0SI * ) ) " f"^ U x U 2 ' + U 2 U x ')3}d 3 A d 3 y. 


In this equation, x has been equated to zero under the argument that the 

separation time x is negligible over the volume of a turbulent eddy. For 

this argument to be reasonable, the convection Mach number should be 

small (see ref. 12, p. 9). The acoustic intensity is seen to consist of 

2 

two primary contributions. The first involves R^ which gives the inten- 
sity which results from the interaction of turbulence with turbulence. 
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This contribution is identified as "self-noise”. The second contribution 
results from the interaction of mean shear (i.e., spatial variation in the 
mean velocity) and turbulence. This contribution is called "shear-noise" 
and can be shown to be non-negligible only in the presence of mean shear. 

The functional behavior of the second-order velocity correlation 
is known from the kinematics of isotropic turbulence (see ref. 10) to be 
the following: 

R xx ■ + "“ir 2 " Tpi < 29 > 

where X = jxj is the distance between point P (where the turbulent 
velocity in the x-direction is u x (y) and point P' (where the turbulent 
velocity in the x-direction is u x (y + X), X x is the component of X_ in 
the x-direction (see fig. 6), and f(X) is a spherically symmetric function 
which vanishes for large values of X. A simple function which has the 
required properties for f(X) is the following ( see ref. 8): 

f(X) = exp (-ttX 2 /L 2 ) (30) 

where L is a characteristic length for the turbulence. With this form for 
f(X), the velocity correlation is given by 

R xx - n [1 - ~ (X 2 - X x 2 )>xp (- —■) . (31) 

L L 

The acoustic intensity in eq. 28 depends on the integral of functions 
of over all values of the separation vector X.' This integration is 

easily performed using spherical coordinates. The self-noise contribution 
to the intensity from a unit volume of turbulence is calculated using 
eq. 31 to be 

2 

8(1 - 4Ncos<J>) R xx 2 (X)d 3 X = L 3 (l - 4Ncos<{>) . (32) 

In order to calculate the shear-noise contribution to the acoustic 
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intensity, the variation of the mean velocity in the jet must be known. For 
weak crossflows the mean velocity is directed primarily in the direction of 
the jet at the exit plane (i.e., in a direction normal to that of the cross- 
flow) for the first eight jet diameters downstream. It is therefore assumed 
that the mean velocity is directed in the x^-direction only. It should be 
noted that this assumption has previously been made in conjunction with the 
eddy convection velocity. In its present context the assumption implies 
the neglect of the mean shear which results from the mean velocity in the 
crossflow direction. It is this mean shear which is responsible for rolling 
the edges of the jet into vortex-like structures. With this assumption, the 
mean velocity in the direction of the observer can be expressed in terms of 
the mean velocity in the x^-direction (i.e., U^) as follows: 

U = ILcose 
x 1 

where 0 is the polar angle between the observer and the x^-axis (see 
fig. 5). This assumption implies the neglect of the terms involving 

and u 2 U 'x in eq " ^ ata Patrick (ref . 13) indicates that 

the mean velocity in the crossflow jet, as in the jet with no crossflow, 
varies rapidly across the jet but slowly along the jet. Ribner's model 
(see ref. 8 ) for the variation of mean velocity in the jet with no cross- 
flow can therefore be used for the crossflow jet. In this model the mean 
velocity is considered to vary rapidly across the jet but no variation is 
included along the jet. The specific form for the mean velocity variation 
is the following: 2 

2 a7rx ? 

u iV = u i ex p (_ * — f-) 

L 

where a is a parameter which is related to the eddy size of the turbulence. 
This approximation is only expected to be valid over the dimension of a 
typical turbulent eddy. The above model for the mean shear gives the follow- 
ing shear-noise contribution to the acoustic intensity from a unit volume of 

turbulence: __ 

TT 2 -2 3 

2 ( 3 IL u L a 4 

(1 - 2Ncos4>) cos 0 R U.U/d A » — r- pr (1 - 2Ncos<j>) cos 0 . (33) 

j xx J. i 2 (1 + a r /Z 
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This contribution to the acoustic intensity is equal to zero when a equals 
zero which corresponds to the case of a uniform flow; therefore, this con- 
tribution is closely associated with the presence of mean shear. 

The results given in eqs. 32 and 33 combine in eq. 28 to give the 
intensity per unit volume of turbulence as follows: 


i(r,e,<j>) /2 2 5 2 


4.-2, %3 
P Q U) f (u ) L 


C(l-4Ncos<)>) + 


U x a 


2 ir a Q r (1-M^cos0)" 


2 3/2 (1+c) 3/2 u 2 


(l-2Ncos(j))cos^03 

where cos0 has replaced x^/r. The factor multiplying the shear- 
noise contribution contains a which is related to the eddy size. As stated 
by Ribner, this factor is difficult to determine but is estimated to be equal 
to unity in the mixing region of the jet. With this estimate, the acoustic 
intensity per unit volume of turbulence is given by 

i(r, 0,<(>) 


4,-2 2 3 
P q to f (u ) L 

0 3/2 2 5 2 

2 v a r 
o 


C(1-4N costfr) + (1-2N cos<j> )cos 
(1 - M^cos0) 3 


(34) 


This equation is the desired result for acoustic intensity in the wind 
tunnel configuration. 


RESULTS 


The acoustic intensity in a moving medium from a unit volume of 
turbulence in a jet flow is given by eq. 34. The medium has a velocity 
a^N in the positive x 2 ~direction , and the convection velocity of a turbu- 
lent eddy is in the positive x^-direction. The directional distribu- 

tion of the acoustic intensity at a constant distance from the jet is given 
by the following expression: 


4 4 

(1 + cos 6) - 2Ncos<t>(2 + cos 0) 

(1 - M^cosO) 3 


(35) 
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where <j> is the polar angle between the wind direction and the observer, 
and 0 is the polar angle between the jet exit direction and the observer 
(see fig. 5). This pattern is plotted on fig. 7 for an observer in the 

7T 

plane of the jet where 4> = y ” The reference intensity is the pattern 

4 z 

given by (1 + cos 0) which is referred to by Ribner as the "basic jet 
directivity". The curve of zero dB therefore indicates the intensity pattern 
which would be emitted by a "jet" with = N = 0. The intensity patterns 
from a jet with M. = 0.54 are presented on fig. 7 for two values of N. 

The pattern given by N = 0 represents the distribution of acoustic inten- 
sity from a straight jet and is symmetric about the x 1 ~axis. The Doppler 
factor (1 - M^cosG) increases the intensity by a large amount in front 
of the jet and decreases the intensity behind the jet. The second pattern 
shown on fig. 7 is the acoustic intensity from a jet in a weak crossflow 
(V./V 2 £ 8). The crossflow Mach number of N = 0.07 produces an asymmetric 
pattern with a greater intensity in the upstream direction and a reduced 
intensity in the downstream direction. Even for this weak crossflow, the 
difference in acoustic intensity measured at corresponding positions upstream 
and downstream is as much as 1.5 dB (at 6 = 315° and = 45°) . As discussed 
previously, this upstream shift results from the convection of signals by the 
acoustic medium. 


Under the assumption of a weak crossflow, the acoustic field which 
is generated by a crossflow jet in motion (i.e., the flight configuration) 
can be calculated directly using the analyses of Lighthill, Ffowcs Williams, 
and Ribner (refs. 5, 6, and 8, respectively). A summary of this calculation 
is presented in the Appendix of this report . In the calculation the velocity 
of the aircraft is a Q N in the negative X 2 ~direction while the turbulent flow 
is characterized by an eddy convection velocity of in the positive 

x^-direction. The directional distribution of the acoustic intensity at a 
constant distance from the aircraft is obtained from eq. A9 to be the 
following: 


V-L > COS UJ 


(1 + Ncos<j>)(l - M^cos 


(36) 


where 0 is the polar angle between the jet axis and the observer position 
and <j> is the polar angle between the positive X 2 ~direction and the observer. 
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The distribution of acoustic intensity for an observer in the plane of the 

7T 

jet (i.e., the plane in which <f> = j- - 0) is presented on fig. 8 . The 

reference intensity is again taken to be the basic jet directivity 
4 

(1 + cos 0). With no aircraft velocity (i.e., N = 0), the intensity pattern 
for the same jet as that on fig. 7 (i.e., = 0.54) is identical with that 

described for the wind tunnel configuration with N = 0. The aircraft Mach 
number of 0.07 shifts the pattern slightly in the direction of motion of the 
aircraft. 

The difference in intensity between the two patterns shown on fig. 7 
is greater than the corresponding difference on fig. 8 . That is, a cross- 
flow of small Mach number in the wind tunnel configuration produces a greater 
change in acoustic intensity than an equal aircraft Mach number in the cor- 
responding flight configuration. This can be demonstrated quantitatively 
when N is small compared to unity (i.e., N<< 1) by expanding the denom- 
inator of eq. 36 to give the following: 

(1 + cos^9) - Ncostj>(T+ cos^ 8 ) ^7 

(1 - M^cosS)^ 


The difference between eq. 37 and eq. 35 is the difference between the 
directional distribution of intensity for the wind tunnel configuration 
and the flight configuration. In the plane of the jet this difference is 

4 

N sin0(3 + cos 8 ) 

_ - — . (.Jo; 

(1 - M^osO) 

This expression is zero for N = 0, 0=0, and 8 = n . For values of 
0 < 8 <ir the intensity in the wind tunnel is less than that for a flight test 
while for 7 r< 0 < 2 ir the intensity is greater for the wind tunnel test. That 
is, the asymmetry of the intensity pattern for a given value of N is greater 
for the wind tunnel configuration than for the flight configuration. This 
results provides a comparison between wind tunnel measurements and measure- 
ments in a flight test for small values of N. 
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CONCLUDING REMARKS 


The result of the analysis for the directivity of the acoustic inten- 
sity from a crossflow jet in the wind tunnel configuration is given by eq. 35. 
This equation can be compared with the corresponding result for the flight 
configuration given by eq. 36. It has been necessary to make several assump- 
tions in deriving these equations which are known to be not strictly valid 
for the crossflow jet. Equations 35 and 36 therefore must be considered to 
provide only a qualitative prediction of the acoustic intensity. The values 
of the constants and parameters which multiply the directional distribution 
of the intensity in eqs. 34 and A9 give at best an order of magnitude estimate 
for the acoustic intensity. It is felt that the directional factor gives a 
reasonable prediction of the directivity from the crossflow jet in the absence 
of reflecting surfaces. This belief in the basic correctness of the directiv- 
ity patterns is based on Ribner's comparisons (see ref. 8, p. 130) of the 
directivity of the straight jet (i.e., no crossflow) with experimental data. 

In this comparison, the patterns of acoustic intensity which are calculated 
using the directional distribution of eq. 35 with N = 0 give reasonable agree- 
ment over most values of 0 with experimental data taken on straight jets. 

For angles near the jet axis, however, the theoretical distribution predicts 
much larger values for acoustic intensity than are observed. Ribner attributes 
this discrepancy to refraction of sound by the mean flow of the jet. In this 
report, refractive effects are not considered by the acoustic theory, and 
similar discrepancies are expected to occur. 

Improvement of this analysis is possible when a more complete picture 
of the crossflow jet is available. One area for improvement is the possible 
refraction of noise emitted by the jet which results from the deviation of 
the crossflow around the jet. Since the jet flow acts as a partial obstacle 
to' the crossflow, there exist gradients in the crossflow velocity near the 
jet flow. This variation in the crossflow may result in a refraction of the 
sound generated in the jet and provide another mechanism for the distortion 
of the sound field. 

It should be recalled that in the calculations of the results, curva- 
ture of the jet by the wind has been neglected. This has been justified by 
restricting the analysis to weak crossflows and postulating that the noisy 
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regions of the jet should be the first eight jet diameters downstream of the 
exit. Curvature of the jet flow, however, could interfere with the propaga- 
tion of signals by shielding some observer positions from signals emitted 
from the jet. 

The results of this analysis for the directional distribution of 
acoustic intensity from a jet in a weak crossflow provide a comparison- for 
corresponding wind tunnel and flight measurements. As demonstrated on ■ 
figs. 7 and 8 and eq. 38, the wind tunnel configuration is more "efficient" 
at shifting the intensity distribution than the flight configuration for a 
given value of N. The effect of the convection of signals by the. acoustic 
medium in the wind tunnel is therefore an important consideration in inter- 
preting acoustic data taken in the wind tunnel tests. 
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APPENDIX 


Acoustic Analysis for the Flight Configuration 


The acoustic analysis presented here is a summary of the contributions 
of Lighthill, Ffowcs Williams, and Ribner. The observer (e.g., a microphone) 
is at rest in a uniform acoustic medium which is also at rest. The source of 
noise is a turbulent flow which can be thought of as being generated by an 
engine of a moving aircraft. This source will be modeled by the flow from a 
cold subsonic jet. 


The analysis of the acoustic field generated by turbulent flows has 
been formally developed by Lighthill (refs. 5, 6). In his development, the 
continuity equation and the momentum equation are combined to give the 
following governing equation: 


a 2 9 a 2 3 2 T.. 

1_£. _ a 2 l.P_ = U 

2 o . 2 3x.3x. 

at ax. i' 2 


(Al) 


where a q is the speed- of sound in the acoustic medium, 
2 

-a p6. . and for a Newtonian fluid, 
o ij 


T. . 


ij 


+ p Vj 


P, . 
ij 


3u. 

■ P5 ij + 11 » sT 
J ] 


3x i 3 9x k 13 


The left-hand side of eq. Al represents the propagation of acoustic waves 
through a stationary-uniform medium. The right-hand side has the form of 
a quadrupole source which is at rest. The quantity T„ is a stress 
produced by pressure, viscous shear stresses, and momentum flux. 


The solution to eq. Al is given as follows: 


p(x,t) - P q = 


, 2 ax.ax. 

4ira o i 2 


|x-y| 

T..(z_,t - 1 — ±L ) 

fiJ a - 


I2L-Z.I 


o ,3 
— d z 


(A2) 
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where is the undisturbed density of the medium, x is the position of 
the observer, and is the position in the turbulent region (see fig. 9). 
The integral is evaluated over the entire region of turbulence. 


The solution given by eq. A2 expresses the density disturbance at 

point x and time t which results from a region of turbulence which is 

at rest relative to the observer. If the turbulence is being generated by 

the engines of an airplane which is in motion relative to the observer, the 

solution must be altered. This situation is considered by Ffowcs Williams 

(ref. 6). A new coordinate system is defined which moves with the aircraft 

at velocity - . In the new coordinate system, a point in the turbulence 

which is moving along with the aircraft is identified by the position vector 

n . The relationship between the coordinate systems is y_ = n_ - < 2 Q Nt when 

both specify the same position in the turbulence. Signals which are received 

I x-y I 

by the observer at time t were emitted at time t - J L . At this time, 

a 

o 


the relationship between the coordinates is X = H - <2 N t + N | x-y] as 

indicated on fig. 9. In terms of the new coordinate system, T. . becomes 

| x~y I 13 

a function of position n_ and the retarded time t — — L . The coordinate 

transformation has the property of changing the volume element in the turbu- 

3 3 

lent region with the relationship between d y. and d n_ given by the 
Jacobian of the coordinate transformation as follows: 



d x Cl + 


N 


(x-y) 


IK-Xl 


. 


The solution for the density disturbance which results from a turbulent 
region of fluid moving with the aircraft is 


p(x,t) - p = 


. 2 3x.3x. 

4ir a x i 

o J 


ij - 


C|x-x| + N.(x-x)] d 11 * 


(A3) 


This solution can be given in a different form by carrying out the differentia- 
tion of the integrand. If the observer position is restricted to the acoustic 
far field of the source, the solution is given by the following: 
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p(x,t) - P Q = - 


4ir a 


(yyptyy.,) 

H|*-y] + N-(x-y)I] :i 3t 


a 2 T. . (n,t ) ,3 
ii — a do 


(A4) 


where t = t - 
a 


| | 

a 


N.(x-y) -.-3 

The Doppler factor LI + — — J represents 


the influence of the movement of the sources of the sound field. 

The expression for the acoustic intensity is obtained from eq. A4 to 
be the following: 


I = 


2 5 

16tt a p 

o o 


( X i - y i ) (x i' y j )(x k~ Z k )(x r Z i ) 

C 1 2_1 + E' + N • (x-z) St^ 2 


3 2 l 1 .(n.t a )3 2 T kt ( i .t 6 ) 


at 


B 


(A5) 


3 3 

do d J ? 


where t = t - 

p 


x-z 


and z^ and have the same relationship as £ and 


n_. The equation for the acoustic intensity can be put into a more useful 

form following the development of Ffowcs Williams. The procedure is identical 

to that discussed previously in this report in conjunction with the wind 

tunnel configuration. A time difference t = t^ - t^ is defined and the 

random turbulence variables are assumed to be stationary functions of time 

in the _n and £ coordinates. A separation vector A_ is introduced which 

is related to r|_ and as follows: • The acoustic intensity is 

then expressed in terms of the correlation function given by S iJH * 

T..(n»t )T, „ (n+A,t +t) as follows: 
xj — a k£ a 


I = 


.,2 5 

16it a p 
o o 


II 


(x -y ) (x -y ) (at -z.) (x^-z^) g 4 33 

2 h 3 —T * 11k « (n,A,x)d JTcTA. 

[|x-zl + N • (x-y.) 3 C | x-zj + N • (x-z) H 9 t ^ 

(A6) 


As pointed out by Ffowcs Williams, the solution in this form is 
not suited for the calculation of the acoustic intensity generated by 
turbulence which is convected relative to the aircraft. The effect of eddy 
convection is included by defining a new separation vector in terms of a 


30 



coordinate system which moves relative to the aircraft at the velocity of 
convection of a turbulent eddy, a M. A new correlation function 
P ijkJ?/— ,T ^ -*- s Refined by expressing A_ in terms of X_ and x with the 

simplification that the separation A_ should be small compared to (x-^) . 

The result of this calculation is given by Ffowcs Williams to be the follow- 
ing equation for the acoustic intensity: 


I = 


1 

i, 2 5 

16ir a p 
o o 


- (x.-y 1 )(x,-y.)(x k -y k )(x r y g ) ^4_ 

' U2L - xl + E‘ (x“Z)^ 5 


P ijkJl ( I I’l.^dWx . 


(A7) 


This equation is simplified further by restricting the observer position to 
distances from the turbulence which are large compared to the dimensions of 
the turbulent flow field. With this restriction, the acoustic intensity is 
given by 


I(r) = 


X.X.X. X. 
1 1 X ft 


16 W (1 + MJS) ( 1 _ MJV 

r r 


3 ,,„(n,x,x) d 3 n d 3 x (A8) 


8x 


4 ijkft — * 


where r = |xj . The acoustic intensity which is calculated from eq. A8 is then 
applicable to the acoustic field generated by the turbulent flow from the 
engines of an aircraft which is traveling at a velocity -q^N through a 
uniform stationary medium. The calculated intensity would be measured by an 
observer at rest in the uniform medium. 


In order to use Ribner's method of modeling the correlation function 
P ijk£^— (ref. 8), the turbulence in the flow from the jet is assumed 
to be convected with a velocity only in the x^-direction. Again, this 

assumption implies that the curvature of the jet flow resulting from the 
crossflow is small (i.e., the crossflow is weak). The velocity of the air- 
craft motion is given by a Q N in the negative X 2 ~direction. With the 
assumption of a weak crossflow, Ribner's model for the correlation function 
is applicable to the crossflow jet. The arguments which lead to this model 
have been outlined previously. The result is the following equation for the 
acoustic intensity which would be emitted by a unit volume of turbulence in 
the mixing region of the jet: 


31 




where <j> is the polar angle between the observer and the X 2 ~axis and 0 is 
the polar angle between the observer and the direction of the convection of 
the turbulent flow (i.e., the x^-axis) . This result provides the acoustic 
intensity per unit volume of turbulence in a jet under flight conditions. 

It provides the comparison to eq. 34 for the intensity from a jet in the 
wind tunnel configuration. 
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FIGURE 1. 


Diagram of a circular jet exhausting into a crossflow. 









FIGURE 7. Directional distribution of acoustic intensity in the plane of 

a i (l+cos^0) - 2 NsinS (2+cos^0) . 

a crossflow jet m a wind tunnel — * E 

(1 - R.cos0) 

(Mj = 0.54; M x = 0.65M^) 
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FIGURE 8. 


Directional distribution of acoustic intensity in the plane of 


a crossflow jet under flight conditions 
(Mj = 0.54; Mj. = 0.65 M ) 


(1 + cos 9) 

(1+N sin6) (1-M^ cos0)' 
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FIGURE 9. 



Diagram showing the relationship between the coordinate 
systems for and n_ at the time when the signal which 

arrives at the observer position x was emitted. 
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